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Abstract 



Power-law sensitivity to initial conditions at the edge of chaos provides a 
natural relation between the scaling properties of the dynamics attractor and 
its degree of nonextensivity as prescribed in the generalized statistics recently 
introduced by one of us (C.T.) and characterized by the entropic index q. We 
show that general scaling arguments imply that 1/(1 — q) = 1/ a m in — 1/ a ma x , 
where a m i n and a max are the extremes of the multifractal singularity spectrum 
f(a) of the attractor. This relation is numerically checked to hold in standard 
one-dimensional dissipative maps. The above result sheds light on a long- 
standing puzzle concerning the relation between the entropic index q and the 

underlying microscopic dynamics. 
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Nonextensivity is inherent to systems where long-range interactions or spatio-temporal 
complexity are present. Long-range forces are found at astrophysical as well as nanometric 
scales. Spatio-temporal complexity, a term introduced to describe the presence of long-range 
spatial and temporal correlations, is found in equilibrium statistical mechanics to emerge 
at critical points for second order phase transitions. Further, the concept of self-organized 
criticality has been recently introduced to describe driven systems which naturally evolve 
through transient non-critical states to a dynamical attractor poised at criticality Long- 
range spatial and temporal correlations are build up in these systems as a consequence of 
avalanche dynamics connecting metastable states. Self-organized criticality is conjectured 
to be in the origin of fractal structures, noise with a 1/f power-spectrum, anomalous dif- 
fusion, Levy flights and punctuated equilibrium behavior 0, which are signatures of the 
nonextensive character of the dynamics attractor. 

The proper statistical treatment of nonextensive systems seems to require a generaliza- 
tion of the Boltzmann-Gibbs-Shannon prescription based in the standard, extensive entropy 
S = —J2iPi^ n Pi ( m units of Boltzmann constant), which provides a link between the mi- 
croscopic dynamics and macroscopic thermodynamics. Inspired in the scaling properties of 
multifractals, one of us |5| has proposed a generalized nonextensive form of entropy 

s q = 1-E f* (I> = i;gerc), (l) 

which recovers the usual entropy form in the limit of q — > 1. The entropic index q controls 
the degree of nonextensivity reflected in the pseudo-additivity entropy rule S q (A + B) — 
S q (A) + S g (B) + (1 — q)S g (A)Sq(B), where A and B are two independent systems in the 
sense that the probabilities of A+B factorize into those of A and of B. A wealth of works has 
been developed in the last few years showing that the above nonextensive thermostatistical 
prescription retains much of the formal structure of the standard theory such as Legendre 
thermodynamic structure, H-theorem, Onsager reciprocity theorem, Kramers and Wannier 
relations and thermodynamic stability, among others [[|. Further, it has been applied to a 
series of nonextensive systems such as stellar polytropes f|, ferrofluids ||, two-dimensional 

2 



plasma turbulence [0], anomalous diffusion and Levy flights ||, cosmology 0, peculiar 
velocities of galaxies ]TIJ and inverse bremsstrahlung in plasma [II], among others ]12| . 



In spite of these variety of applications of nonextensive thermostatistics, a full and gen- 
eral understanding on the precise relation between the entropic index q and the underlying 
microscopic dynamics was still lacking. It has been conjectured that the generalized ther- 



mostatistics is a natural frame for studying fractally structured systems |12j and simple 
relations were found between q and the characteristic exponents of anomalous diffusion and 
Levy flights distributions [§]. Further, recent works have shown that the entropic index q has 
a monotonous dependence on the fractal dimension df of the dynamical chaotic attractor of 



dissipative nonlinear systems [13 



The aim of the present work is to develop the precise connection between the nonexten- 
sivity parameter q and the scaling properties of the critical attractor of nonlinear dynamical 
systems. Particularly, a prototype complex dynamical state will be taken to be the onset 
of chaos of nonlinear low- dimensional maps as it is well stablished that this state displays 
long-range temporal and spatial correlations. We will show that the power-law sensitivity 
to initial conditions at the edge of chaos provides a natural link between the entropic index 
q and the attractor's multifractal singularity spectrum. By using typical scaling properties 
of Feigenbaum-like maps we will show that 1/(1 — q) = l/a m i n — l/a max , where a m i n and 
&max are the scaling exponents related respectively to the most concentrated and most rar- 
efied portions of the attractor. We numerically illustrate the above result using standard 
one-dimensional maps such as generalized logistic-like maps and the circle map. 

For the sake of simplicity, we will concentrate our attention to the simple case of one- 
dimensional non-linear dynamical systems. One of their most prominent features is related 
to the sensitivity to initial conditions. In order to quantify this aspect, Kolmogorov and 
Sinai's definition of the rate at which the amount of information about the initial conditions 
varies can be seen as : 

K x = lim Jim lim ^-[S^N) - S l (0)}, (2) 
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where r is a characteristic time step (in fact, r — > for differential equations; r = 1 for 
discrete maps) and Si(0) and Si(N) stand, respectively, for the entropies evaluated at the 
times t — and t = Nr. The entropy can be evaluated by considering an ensemble of 
identical copies of the system and defining pi as the fractional number of copies that are 
in the i-th cell (of size I) of the phase space. If one uses the extensive Boltzmann-Gibbs- 
Shannon entropy form S = Si = — X^iPjhipj, where W is the number of configurations 
at time t, the Kolmogorov- Sinai entropy results in 

K x = limlim lim —lnW(N)/W(Q), (3) 

where equiprobability (pi = 1/W) was assumed. Notice that the above expression implies 
in an exponential sensitivity to initial conditions W(N) = W(Q)e 1 . K\ plays (consis- 
tently with Pesin's equality) the role of the Liapunov exponent Ai which characterizes the 
exponential deviation of two initially nearby paths = limAa;(o)->o Ax(t)/Ax(Q) = e Xlt 
(£(£) is solution of d£/dt = Ai£). When Ai < (Ai > 0) the system is said to be strongly 
insensitive (strongly sensitive) to the initial conditions. The marginal case of Ai = occurs 
at the period doubling and tangent bifurcation points, as well as at the cumulating point 
of the period doubling bifurcations, i.e., at the threshold to chaos. The failure of the above 
scheme in distinguishing the sensitivity to initial conditions at these special points is related 
to the nonextensive (fractal-like) structure of their dynamical attractors. 

Recently it was argued that, within the generalized nonextensive entropy of Eq. 1, the 



sensitivity to initial conditions of one-dimensional nonlinear maps becomes expressed as JT3 

e(*) = [l + (l-ff)V] 1/(1 ~ a) . ( 4 ) 

which is solution of d£/dt = \ q !; q . The above expression recovers the usual exponential 
sensitivity in the limite of q — > 1 (extensive statistics). Further, it implies a power-law 
sensitivity when nonextensivity takes place (q ^ 1). Previous numerical calculations have 
shown that the period doubling and tangent bifurcations exhibit weak insensitivity (q > 1) 
to initial conditions (T^j. At the onset of chaos, weak sensitivity (q < 1) shows up and 



the value of q was numerically verified to be closely related to the fractal dimension df of 
the dynamical attractor [[Uj. In what follows we will use scaling arguments to analytically 
express the entropic index q as a function of the fractal scaling properties of the attractor. 

Actually, as in many nonlinear problems, the scaling behavior of the attractor is richer 
and more complex than is the case in usual critical phenomena. It is necessary to introduce 



a multifractal formalism in order to reveal its complete scaling behavior ||14|| . A central 
quantity in this formalism is the partition function Xq{N) = J2iLiPi, where pi represents 
the probability (integrated measure) on the i-th box among the N boxes of the measure 
(we use q instead to the standard notation q in order to avoid confusion with the entropic 
index q). For example, in chaotic systems pi is the fraction of times the trajectory visits 
the box i. In the iV — > oo limit, the contribution to Xq{N) oc iV~ r ®, with a given q, 
comes from a subset of all possible boxes, whose number scales with N as Nq oc N^^\ 
where f(q) is the fractal dimension of the subset (f(q = 0) is the fractal dimension df of 
the support of the measure). The content on each contributing box is roughly constant 
and scales as Pg oc N~ a ^\ These exponents are all related by a Legendre transformation 
r (<?) = q a {(i)~ /(?)■ The multifractal object is then characterized by the continuous function 
f(a), which reflects the different dimensions of the subsets with singularity strength a. 

The multifractal formalism has been widely used to characterize the spectrum of singu- 
larities of some important objects arising in nonlinear dynamical systems as for example the 
critical cycle at the onset of chaos |L4| , diffusion limited aggregation, fully developed turbu- 



lence, random resistor networks, among many others ||15|| . The a values at the end points of 
the f(oe) curve are the singularity strength associated with the regions in the set where the 
measure is most concentrated (a m i n = a(q = +00)) and most rarefied (a max = a(q = —00)). 

The scaling properties of the most rarefied and most concentrated regions of multifractal 
dynamical attractors can be used to estimate the power-law divergence of nearby orbits. 
Consider the set of points in the attractor generated after a large number B of time steps 
(Pi = 1/B is therefore the measure contained in each box). The most concentrated and 
most rarefied regions in the attractor are partitioned respectively in boxes of typical sizes 



/ +00 and l-oo- From these, one may determine the end points of the singularity spectrum 
as a min = \npi/lnl +00 (hence l +OQ oc B~ 1/ami ") and a max = lnpi/lnLo, (hence Z_oo oc 
£>-i/a max y Further, the smallest splitting between two nearby orbits, which is of the order 
of l +OQ , becomes at most a splitting of the order of l-^. With these scaling relations Eq. 4 
reads: 

Loo/Z+oo oc (5) 
which implies a precise relation between the entropic index q and the extremes of f(a): 

1 1 1 (6) 



1 Q O-min C^max 

The above expression is the main result of the present work. It asserts that, once the 
scaling properties of the dynamical attractor are known, one can precisely infer about the 
proper nonextensive statistics that must be used in order to predict the thermodynamics of 
the system. 

Let us illustrate the above result using as prototype multifractal objects the critical 
attractor of one-dimensional dissipative maps. The scaling properties near the onset of 
chaos allows for analytical expressions for the end points of the singularity spectrum. As it 
has been shown by Feigenbaum, with B = 6 n -cycle elements on the attractor (b stands for 
a natural scale for the partitions), the most rarefied and most concentrated elements scale, 
respectively as l-^ ~ a~p n and l +00 ~ a~p zn , where a? is the Feigenbaum's universal scaling 
factor and z represents the nonlinearity (inflexion) of the map at the vicinity of its extremal 
point fl6|| . Since the measures there are simply p_oo = p +oa = Pi = b~ n , these end points 



are respectively 

- lnb (7\ 
map 

In b 

Therefore, the entropic index q can be put as a function of the Feigenbaum's scaling factor 
as 



In order to determine the singularity spectrum of the critical attractor, we implemented 
the algorithm proposed by Halsey el al |14]], which has been shown to be more precise than 



the standard box counting method. Firstly, we consider a family of generalized logistic maps 

x w = 1 - a\x t \ z (1 < z < oo; < a < 2; - 1 < x t < 1), (10) 

which exhibits a period-doubling cascade accumulating at a c (z) (b = 2 is therefore the 
natural scale for the partitions). Here z is precisely the inflexion of the map at the vicinity 
of its extremal x — 0. Typical multifractal singularity spectra are shown in Fig. 1. From 
their end points we can estimate the z-dependence of the universal scaling factor a.p(z). 
The so obtained values are shown in figure 2 together with known asymptotics fl7j and 
the parametric dependence of a m i n and a max on the fractal dimension df (see inset). The 
entropic index q was independently obtained (numerically) from the plots of ln£(iV) = 
J2tLi l n [a^l^t Y^ 1 } versus IniV, where N is the number of iterations. The upper bound of 
these plots have slopes equals to 1/(1 — q) (see Fig. 3a). Its fractal- like structure reflects the 
presence of long-range temporal correlations at the critical point. The so obtained values of 
q are plotted in Fig. 4 against the numerical values l/a m i n — l/a max and corroborate the 
relation predicted from scaling arguments. 

We also computed the multifractal spectrum f(a) and the sensitivity function £(N) for 
the following two-parameter map: 



x t . 



^ = dcos(n\x t - l/2\ z/2 ) (1 < z < oo; < d < oo; -d<x t <d). (11) 



This map also display a period doubling route to chaos (z is the inflexion of the map at the 
vicinity of the extremal point x = 1/2). A typical onset to chaos is found, for example, to 
be at d c {z = 2) = 0.865579... . Our numerical results confirm, as expected, that this map 
has, for fixed z, the same f(a) and q of the logistic-like map, i.e., both maps belong to the 
same universality class. 



As a final illustration, we computed £(iV) for the circle map which is an iterative map- 
ping of one point on a circle to another. This map also exhibits a transition to chaos via 
quasiperiodic trajectories. It describes dynamical systems possessing a natural frequency uj\ 
and driven by an external frequency LO2 (fi = 0J1/0J2 is the called bare winding number) and 
belongs to the same universality class of the forced Rayleigh-Benard convection ||TB| . These 
systems tend to mode-lock at a frequency uj\ (the ratio between the response frequency to 
the driving frequency uo* = uj\juj2 is usually called dressed winding number). The standard 
one-dimensional version of the circle map reads: 

e t+l = t + n- — sin(27T0 t ), (0 < Q < 1; < K < 00; < 6 t < 1). (12) 
2n 

This map exhibits a nonlinear inflexion around its extremal point only for K > 1. Therefore 
K = 1 is the onset value above which chaotic orbits exist (for K < 1 the orbits are always 
periodic). A well-studied transition takes place at K = 1 and dressed winding number equal 
to the golden mean u 

* = <y<5- l)/2 = 0.61803..., which corresponds to Q = 0.606661... . 
With these parameters, the map has a cubic inflexion [z = 3) near its extremal point 
9 = and its universal scaling factor is found to be ap = 1.2885... [[L9| (b = is the 
natural scale for the partitions). From Eq. 9 the predicted value for the entropic index is 
q = 0.0507.... This value also agrees with the numerical estimation based on the sensitivity 
to initial conditions (see Fig. 3b). 

In conclusion, we have shown how the proper nonextensive statistics can be inferred from 
the scaling properties of the dynamical attractor at the onset of chaos in one-dimensional 
dissipative maps. The relation between the entropic index q of generalized statistics and the 
multifractal singularity spectrum of the dynamical attractor was derived using quite general 
scaling arguments applied to the most concentrated and rarefied regions of the attractor. The 
proposed relation is therefore expected to hold for higher dimensional dissipative systems 
and to provide a close relationship between the nonextensive statistics formalism and the 
self-organized critical states of large driven dynamical systems. Analogous connections might 
exist for Hamiltonian systems with long-range interactions. 
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I. FIGURE CAPTIONS 



Figure 1 - Multifractal singularity spectra of the critical attractor of generalized logistic 
maps with z = 1.5,2.0 and 3.0 as numerically obtained following the prescription in Ref. 

0- 

Figure 2 - The logistic-like- map values of ap{z) (numerically obtained from both a m i n 
and a max ). The solid lines represent known analytical expressions for the asymptotic be- 
haviors ([ap(z)] z — > 1/0.033381... as z — > oo; and chf(z) ~ — l/[(z — 1) In (z — 1)] as z — > 1 
|T7|]). [ Inset: a min and ot max versus df. Dashed lines are to guide the eyes. 

Figure 3 - ln£(iV) versus IniV. (a): standard logistic map (z=2). The solid line 
represents the theoretically predicted slope for the upper bounds 1/(1 — q) = In of(2) / In 2 = 
1.3236... (a F {2) = 2.5029... |l|), hence q = 0.2445.... (b): circle map at K = 1 and 
to* = (v5— 1)/2 (hence f2 = 0.606661...). The solid line represents the theoretically predicted 
slope for the upper bounds 1/(1 — q) = 2 hxatp/ lnu;* = 1.0534... (op = 1.2885... pJ[|), hence 
g = 0.0507.... 

Figure 4 - g versus l/a m i n — l/a max for the generalized logistic map (circles) and for the 
circle map (square). The straight line represents the scaling prediction. 
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